Quantum impurity coupled to Majorana edge fermions 
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We study a quantum impurity coupled to the edge states of a two-dimensional helical topological supercon- 
ductor, i.e., to a pair of counter propagating Majorana fermion edge channels with opposite spin polarizations. 
For an impurity described by the Anderson impurity model, we show that the problem maps onto a variant of the 
interacting resonant two-level model which, in turn, maps onto the ferromagnetic Kondo model. Both magnetic 
and non-magnetic impurities are considered. For magnetic impurities, an analysis relying on bosonization and 
the numerical renormalization group shows that the system flows to a fixed point characterized by a residual In 2 
entropy and anisotropic static and dynamical impurity magnetic susceptibilities. For non-magnetic impurities, 
the system flows instead to a fixed point with no residual entropy and we find diamagnetic impurity response 
at low temperatures. We comment on the Schrieffer- Wolff transformation for problems with non-standard con- 
duction band continua and on the differences which arise when we describe the impurities by either Anderson 
or Kondo impurity models. 

PACS numbers: 73.20.-r, 72.10.Fk, 72.15.Qm, 75.30.Hx 
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I. INTRODUCTION 



Majorana fermionic operators are characterized by the re- 
lation rp = ?/, i.e., they represent particles which are their 
own antiparticles. A pair of Majorana (real) operators can be 
combined into one standard Dirac (complex) fermionic cre- 
ation and one corresponding annihilation operator which op- 
erate within the Fock space of a single-level system. Majo- 
rana fermions can also be thought of as an equal linear com- 
bination of particles and holes. In recent years, a number 
of proposals have been advanced for physical realizations of 
condensed-matter systems having low-energy excitation spec- 
tra which can be formally described using Majorana fermion 
operators (see Refs. QHsHa for recent reviews). In such in- 
teracting many-particle systems Majorana fermions appear at 
low temperatures as emergent degrees of freedom. In particu- 
lar, it has been argued that topological superconducting phases 
can be induced in topological insulators^— by the proximity 
effec t 4 ' 5 ' 12 ^ . These phases are characterized by the presence 
of one or several Majorana excitation branches which are lo- 
calized along the perimeter of the two-dimensional supercon- 
ducting sheet. Such Majorana bands have energies inside the 
superconducting energy gap. In the absence of defects they 
are, in fact, the only excitations inside the gap. One distin- 
guishes chiral topological superconductors with a single Ma- 
jorana continuum which propagates in a unique direction (chi- 
rality) and helical topological superconductors with two Ma- 
jorana continua which counter propagate and have opposite 
spin polarization (helicity or spin-momentum locking) 4 *!^^. 
While there is presently no confirmed physical realization of a 
helical topological superconductor, several recent theoretical 
proposals are based on heterostructures made from topologi- 
cal insulators (such as Bi2Se3) and conventional superconduc- 
tors (such as elemental superconducting metals). By tuning 
magnetic doping of the topological insulating layer, it should 
be in principle possible to obtain both chiral and helical topo- 
logical superconductors^. 



An important issue in this field is the fate of the Majorana 
fermions when electron-electron interactions are taken into 
accoun t 22 ' 23 and the role of disorder-2 4 ^— , in particular that of 
impurities 2 ^. Dilute concentrations of impurities would lead 
to various measurable low-temperature anomalies due to the 
Kondo effect, while higher concentrations could even destabi- 
lize the edge states and qualitatively affect the behavior of the 
system. It is thus of considerable interest to accurately ana- 
lyze the behavior of interacting impurities embedded near the 
edges of topological superconductors (TSC). This is particu- 
larly important because the doping of the topological insulator 
layer by magnetic impurities has been proposed as one of the 
methods to tune their material properties and to control their 
ground states. It is therefore likely that some magnetic dopant 
atoms will be invariably present at the edge area of the TSC, 
where the Majorana states are localized. 

Since the impurities hybridize with the conducting states 
which have half the degrees of freedom of Dirac fermions 
(electrons), it is expected that there will be a number of partic- 
ularities compared to the behavior of standard impurity mod- 
els, such as the Anderson and Kondo impurity models. Fur- 
thermore, it is conceivable that one could attach quantum dots 
to such materials so that they would hybridize with the edge 
states of TSCs; such hybrid structures can also be described 
using the same class of impurity models. 

In this work we study the Anderson impurity model with 
modified conduction-band and hybridization terms which 
only involve half the degrees of freedom of the standard 
model. In Sec. [II] we describe the model and perform map- 
pings to other known models. In Sec. [TTT1 we use bosoniza- 
tion techniques in order to gain some qualitative understand- 
ing about the expected behavior of the model. In Sec. [IVj 
this analysis is complemented by a numerical analysis using a 
reliable non-perturbative technique, the numerical renormal- 
ization group (NRG), which allows to address both thermo- 
dynamic and dynamic properties of the system. By help of 
these two complementary studies, we show that the model ad- 
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mits two qualitatively different types of low-temperature fixed 
points, one associated with magnetic impurities (those with 
large electron-electron repulsion and occupancy near half- 
filling) which is characterized by In 2 residual entropy and 
anisotropic magnetic response, and another associated with 
non-magnetic impurities (those with small electron-electron 
repulsion and/or occupancy away from half-filling) which is 
characterized by zero residual entropy and diamagnetic re- 
sponse. This is different from the standard Anderson impurity 
model where the low-temperature fixed point always belongs 
to the same family of Fermi-liquid fixed points and there are 
no discontinuities as a function of the model parameters (in- 
teraction strength, level energy, hybridization strength), only 
the quasiparticle scattering phase shifts are smoothly changed. 
In Sec.[VJ we perform the Schrieffer- Wolff transformation in 
the case of Majorana bands. We comment on the relations be- 
tween the Anderson and Kondo impurity models in this gener- 
alized setting and compare our results of a magnetic impurity 
coupled to Majorana edge states with a recent study^ 7 which 
uses instead a Kondo Hamiltonian description. In particular, 
we find different results for the dynamical magnetic response 
and observe that the phase transition predicted in Ref. [271 does 
not show up when the impurity is instead properly described 
by a Anderson model. This is due to the fact that the phase 
transition predicted in Ref. |27l occurs for parameter values of 
a Kondo model which cannot arise from the Schrieffer- Wolff 
transformation of an Anderson model. 



II. MODEL AND MAPPINGS 




Figure 1. (Color online) Schematic representation of the single im- 
purity problem studied in this work. A single impurity level (cir- 
cle) is hybridized with two continua of counter-propagating Majo- 
rana fermions with opposite spin polarizations (dashed lines). 

We study the single-impurity Anderson model (SIAM) con- 
sisting of an impurity (single orbital) coupled by hybridization 
to the edge states of a helical topological superconductor, i.e., 
to two continua composed of counter propagating Majorana 
particles with opposite spin polarization. The spin quantiza- 
tion axis is taken to be along the z-axis; the rotation symmetry 
in the spin space is broken in this model. The Hamiltonian can 
be obtained from the standard SIAM by projecting out half of 
the conduction-band degrees of freedom. A schematic repre- 
sentation is shown in Fig. Q] The Hamiltonian is 



where Hi is the impurity Hamiltonian 

Hi = ^(n, - 1/2) + U(n t - l/2)( % - 1/2) 

a 

+B x (l/2)(dl^ + d\d l ) 

+ By(l/2)(id\d t -id\d i ) 

+B z (l/2)(n t -n i ), 



i?2 is the coupling Hamiltonian 



H 2 = Y,V k 



r\ko 



(2) 
(3) 



(4) 



and i?3 is the conduction-band Hamiltonian for the Majorana 
modes 

#3 = ( vkl l-k,tVk,t + v(-k)r)-k,iVk,i.) ■ (5) 

fc>0 

The operators d\ with a I create an electron at the impu- 
rity level, and n a = d\d a is the electron number (occupancy) 
operator. The parameter 5 = e + U/2 measures the departure 
from the particle-hole symmetric point, e is the on-site energy, 
U is the electron-electron repulsion parameter, while B x , B y , 
and B z are the components of the external magnetic field ap- 
plied on the impurity. The r\^,a are Majorana fermions opera- 
tors satisfying 77 j. ff = ?7_ fc:CT and {r]-k, a , Vk',<r'} = 4fc"W- 
v is the Fermi velocity. The hybridization is given by the ma- 
trix elements Vk- As common in the treatment of impurity 
models, we neglect the /s-dependence of Vk, i.e., Vk = V; 
this is permissible since we are mostly interested in the low- 
energy behavior of the system. We furthermore assume that 
Vk is real (see AppendixlAlfor a generalization to complex Vk 
where we show that the phase of Vk determines the orientation 
of the privileged axis, to be discussed below, in the xy plane 
that is perpendicular to the spin-quantization axis of Majorana 
bands, i.e., the z-axis). The hybridization of the impurity with 
the conduction band is fully described by the hybridization en- 
ergy scale r = irp\V\ 2 , where p is the density of states in the 
conduction band. In this work, we consider a constant den- 
sity of states p = 1 Jixv so that T is a constant independent of 
energy. 

It is important to note that in writing the Hamiltonian in the 
form of Eqs. (Q][5]l we have omitted terms which describe the 
coupling of the impurity to the bulk states of the supercon- 
ductor which have energies above the superconducting gap. 
These modes can renormalize the effective parameters of the 
impurity Hamiltonian. For simplicity, we assume that these 
high-energy modes have already been integrated out and that 
the Hamiltonian H is an effective low-energy Hamiltonian 
valid on energy scales below the superconducting gap. The 
bandwidth of the Majorana modes thus equals twice the su- 
perconducting gap. 

We now introduce the local Majorana operators r)i a through 



H = H 1 +H 2 + H 3 



(1) 



Via = -j={d a + dl) , 7]2a = ^ (da - 4) , 



(6) 
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(note that we follow the normalization convention i]f a = 1/2). 
The local spin and isospin operators can be expressed as 

s x = 2 ( rf I d t + d \ d i) = \ (VilV2f + J71TT24-) 
s z = - (d\df - 4 d i) = \ (-VUV21 + VitV2f) 

1 / \ 7 

= 2 (44 + d ± d y = 2 ^ u??2t ~ m tV2i) 

h = 2 (*44 + *^ d t) = | faufif - %4.%t) 

i* = - (d\d t + d\d x - lj = ^ (7714.7724. + 77if?72t) 

We rewrite if 1 as 

ifi = <^y^ ir\\cjf]2a + Uintvwn^i 

<7 

+ B x (»/2) (77 U 772t + >7it^) (8) 

+ B z (i/2) (mtmt - mimi) , 

while 

H 2 =J2vV2tVkaV2a- (9) 

fecr 

Only the two r\2a Majorana local modes are coupled to the 
continuum (hybridized). For 6 = B = U = 0, the two 7710- 
modes are fully decoupled from the rest of the system. For 
U ^ 0, the rjicr modes interact via the quartic term with the 
7720- modes. As we show in the following, this brings about 
non-trivial effects. For 5 = B = and finite U, we will see 
that on the temperature scale T\ = f(U,T) the fluctuations of 
the 772,7 modes are frozen out and that the entropy is reduced 
by 2 x (1/2) In 2 = In 2. The modes remain active, and 
the system has In 2 residual entropy down to T = 0. We 
would therefore like to study the nature and the dynamics of 
this residual degree of freedom. We find that 

*'7H ? /it = s v + V ( 10 ) 

The residual degree of freedom thus corresponds to a mixed 
spin-isospin mode associated with a linear combination of 
spin and isospin projection operators along the y direction (see 
Appendix [A] for a generalization to an arbitrary phase factor 
in Vfc). For this reason we expect that the model has different 
magnetic response in the y direction as compared to the x and 
z directions. In the following, the privileged y direction will 
be referred to as the "longitudinal" direction, while the other 
two will be called "transverse" directions. 

We now map the problem onto an interacting resonant 
two-level model. For that purpose, we introduce new Dirac 



fermionic operators defined by 



afc = 


1 


mi), 


4 = 


1 

71 {71 - kt 


- iv-n) 


-ib = 




imi) , 


1 = 


71 (r/lt + 


ivn) ■ 



(11) 



In terms of these operators, 

Hi =5 (6 f / + H.c.) - U(n h - l/2)(n f - 1/2) 

+B x /2 (iftt/t + H.c.) + B z /2 (&t/t + Re.) , 
+B y /2(l-n b -n f ) (12) 

H 2 = VV2 (4,6 + H.c.) , (13) 

k 

H 3 = ^2vkala k . (14) 

k 

Notice that the summation in the last equation is indeed for all 
k. Here nt = b^b and n/ = f'f are the occupancy operators 
for orbitals b and /, respectively. This is a variant of the inter- 
acting resonant two-level model for spinless Dirac fermions 
with attractive charge-charge coupling. The b mode directly 
hybridizes with the continuum, while the / mode is fully de- 
coupled for S = U = B = 0. We note the similarity between 
the B x and B z terms in Eq. (fT2l) (they differ essentially only 
by a phase shift of the operators), while B y has a qualitatively 
different form. This is in line with the expected anisotropy 
of the system. The sum of the terms in B x and B z can be 
compactly rewritten as (B± b'f' + H.c.)/2, where 

B ± =B z +iB x . (15) 

The interesting behavior results for non-zero parameters S, U, 
and/or B. For finite interaction U, we expect to observe or- 
thogonality catastrophe physics due to the effect of the occu- 
pancy changes at the impurity level / on the continuum. 

In order to capture the physics associated with this model, 
let us first use a field-theoretical approach. This is the subject 
of the next section. 



m. BOSONIZATION APPROACH 

A Hamiltonian somewhat similar to Eqs. (1 1 2H 14b has been 
studied recently in the context of the charging of a narrow 
level in a quantum dot capacitively coupled to a broader 
one^^i. However, contrary to the Hamiltonian in Eq. ( fT2b . 
the levels in the dot were interacting via a repulsive Coulomb 
interaction, in which case the system always has a non- 
degenerate ground state and is a conventional Fermi liquid. 
Here, the interaction is attractive which, as we will see, can 
significantly change the physics. 

We start with the case B y = 0. We follow the bosonization 
procedure developed in Ref. [3^ and repeat the main steps for 
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completeness. We first diagonalize the Hamiltonian H2 + H3 
which corresponds to a resonant-level model. The density of 
states of the level b corresponds to a Lorentzian of finite width 
r centered around zero energy. If we assume 6, \B±\ <C I\ 
we can replace the density of states of the level b by a flat 
density 30 for energies u> € [— ttT/2, ttT/2] and elsewhere, 
with height 1/(ttT). In this limit, the energy scale V plays 
the role of an effective bandwidth. Therefore, H reads in the 
continuum and low-energy limit 

00 

H = -iv J ip ] \x)d x il>{x)dx - Ua : ^(0)^(0) : (n f - 1/2) 

— OO 

+S^(^(0)f + H.c.) + B ± ^E(^(0)p + H.c), (16) 

where a ~ 2v/T is a new UV cut-off associated with the 
new effective bandwidth T, ip(x) is a right-moving field as- 
sociated with the diagonalization of the low-energy limit of 
H-2 + H3. The Hamiltonian in Eq. ( TTol ) can be treated using 
bosonization. We introduce a bosonic field associated to the 



chiral fermionic one ip(x) 



2-Ka 



-i<t>(.x) such that 



f ^( V 0) 2 -i;^V0(O)(n / -l/2) (17) 

4.7T 7T 



where 



ay = arctan([//2r) 



(18) 



is fixed so that the phase shift we obtain for S = B± = in 
the fermionic representation equals the one calculated in the 
bosonic representation for a given occupation of the / level^ 
Next, we apply the canonical transformation H' = WHU 
with 



[/ = exp[i(2aa/7r)0(O)(n/-l/2)], 



(19) 



such that 



00 

H' = v J g(V0) 2 + -^(e^(°)/ + ff.c.) 



+ £X( e <7-*(0) / t + J T>c .) j 

where we have introduced 

7± = 1 ± 2a[//7r. 



(20) 



(21) 



Let us first consider the following two situations: 8 7^ 
0, B± =0 and 5 = 0, B± ^ 0. In both cases, the Hamil- 
tonian in Eq. d20l i maps to a bosonized formulation of the 
anisotropic Kondo model 



where = vk and cr % are the Pauli matrices (i = x, y, z). 
The Kondo couplings Ji are such that J x = _X jj = J±. The 
bosonization of d22L as derived in Refs. 32ll33l provides the 
identification 



7i 



= J±/V8 

1— 2 
= \/2[l arctan(7rpJz/4)] 



(23) 



where p is the density of states in the large-T limit we are 
considering. In the limit where V is much larger than U, we 
can linearize the arctan functions and estimate 



J z w 7iT(2 - V2) =F V2U. 



(24) 



The upper sign ( — ) corresponds to the case (5^0, B± = 0, 
and the lower sign (+) to the case (5 = 0, B± 7^ 0. Let us now 
discuss both cases separately. 

Case 5^0 and B± = 0. - We see from Eq. d24ji that a 
large U ^> T enforces J z < 0, i.e., a ferromagnetic Kondo 
model, while a small U leads to an antiferromagnetic Kondo 



model. In the latter case, we can extract a Kondo scale T 



K 



associated with the screening of the magnetic moment. From 
the known results on the anisotropic Kondo model, 34 we find 



± K 



(25) 



with 7+ defined in Eq. (f2Tb . We remind that this result makes 
sense only for J z > (small-[/ case). 

In order to estimate the low energy physics in the presence 
of both <5 and U terms, we can rely on the lowest order renor- 
malization group (RG) equations for the Kondo model 33 

^±=pJ±Jz-p 2 J±(Jl + J')/4 (26) 

^= p Jl-p2 Jz Jl/2, (27) 

where / is a running length scale. In particular, if we ne- 
glect the cubic terms, this is the standard Kosterlitz-Thouless 
RG flow where J z — J\ is a constant of the flow. When 
J z < 0, the infrared fixed point depends on whether | J z (a) | is 
above or below the separatrix J± = —J z , assuming J± > 0. 
If |J 2 (a)| > J±(a), the fixed point is a strong-coupling 
one which corresponds to the antiferromagnetic Kondo fixed 
point and the formation of a Kondo singlet. However, for 
|Jz(a)| < Jj_(a), the flow is driven toward a fixed line cor- 
responding to J± — >• and J z — > J* < associated with 
the ferromagnetic Kondo model and therefore an unscreened 
free moment. Although we have assumed T to be larger than 
other scales, the condition J± (a) > \ J z (a) \ associated with 
a Fermi-liquid fixed point and a screened magnetic moment 
translates into 



6>U/2 



r. 



(28) 



K 



E 

ha 



E 



(22) 



Case 6 — and B± 7^ 0. - This case is simpler since 
J z > for all values of U . Therefore, the low energy physics 
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is described by a Fermi liquid fixed point with a screened mag- 
netic moment and we can define a Kondo scale TZ by 

fB x \^7E 
T k~ t [-y) - ( 29 ) 

with 7_ defined in Eq. ( f2TT >. 

General case. - When both <5 and B± are different from 
zero, the mapping to the anisotropic Kondo problem breaks 
down. However, at small U compared to T, both terms inde- 
pendently favor a strong-coupling fixed point with a screening 
of the local moment. Therefore, we expect the same behav- 
ior in the presence of both terms. In the large-[/ limit with 
B± = 0, a ferromagnetic Kondo model is favored and a lo- 
cal moment is stabilized. The addition of a small magnetic 
field B± around this fixed point is, however, a relevant per- 
turbation which destabilizes this behavior and a non-magnetic 
fixed point is also expected. 

In our analysis, we have considered B y = 0. However, 
we can extend the previous bosonization analysis to models 
with By 0. The term B y has two effects: First, it provides 
an additional potential scattering term —B y a/2 : ip(0)ip(0) : 
in the continuum-limit Hamiltonian in Eq. ([Tol l which will 
modify the phase-shift analysis. Second, it adds an extra term 
—B y (rif — l/2)/2 which corresponds to a magnetic field of 
the form —B y S z /2 in the equivalent Kondo Hamiltonian in 
Eq. (l22l . As discussed above, the latter term would polarize 
a free local moment and no residual entropy is expected at 
energy below this magnetic field. 



IV. NUMERICAL RENORMALIZATION GROUP 
ANALYSIS 

We now study the Hamiltonian in its original (non- 
transformed) form of Eq. (Q]i using the numerical renormal- 
ization group (NRG) method^—. This technique is essen- 
tially an iterative exact diagonalization of an appropriately 
discretized form of the Hamiltonian, where at each partial di- 
agonalization step the high-energy excitations are discarded 
("truncated") and only the low-energy part of the spectrum 
is retained in the subsequent calculation steps. This approx- 
imation is appropriate because the matrix elements coupling 
levels with widely different energy scales are small^. The 
NRG method has already been applied to problems involving 
an impurity coupled to Majorana edge fermions; for details, 
see Ref. [38l (In this work, the NRG calculations have mostly 
been performed using parameters A = 3, -Bcutoff = 10 or at 
most 5000 states, whichever is lower, N z = 4, a = 0.6.) 



A. Thermodynamic properties: particle-hole symmetric case 

We first study the problem for a particular choice of 5 = 0, 
i.e., we focus on the particle-hole (p-h) symmetric problem. 
We will show in the following that there is actually a finite re- 
gion in the (5, U, T) parameter space where the problem flows 
to the same type of the low-temperature fixed point as along 
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Figure 2. (Color online) Thermodynamic properties (impurity en- 
tropy and impurity susceptibilities) at the particle-hole symmetric 
point S — for a range of values of the interaction strength U. The 
arrows indicate the direction of increasing U. In the absence of the 
external magnetic field, Olie ha.S \xx,imp — X,zz,imp due to symme- 
try. The out-of-diagonal magnetic susceptibilities are all zero. On 
the curves for U = 0.02, we label the positions of the temperatures 
T\ and T2 (see the main text for their definitions). 

the p-h symmetric line, thus the results are more general. The 
impurity entropy and thermal spin and charge susceptibilities 
are shown in Fig.[2]for a range of the interaction strengths U, 
including for the non-interacting U = limit. We remind the 
reader that the impurity entropy quantifies the number of the 
effective degrees of freedom p on the impurity via 

S imp (T) = k B In p(T). (30) 

The impurity magnetic susceptibility corresponds to the effec- 
tive local magnetic moment as 

tf(T) = T X u, imp (T), (31) 

where 

Txij,im P (T) = (SiSj) - (SiSj)o, (32) 

with Si the i-component of the total spin of the system. Here 
()o denotes the expectation value for a system without the im- 
purity. Finally, the impurity charge susceptibility x<2,imp is 
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defined through 



TXQM P (T) = (Q 2 ) - (Q 2 ) 



(33) 



where Q is the total charge in the system. Since the impurity 
susceptibility is defined as the difference between an expecta- 
tion value in the full system and the expectation value in the 
system without the impurity (i.e., the impurity contribution to 
the total system susceptibility), the impurity susceptibilities 
can be negative. This is different from the local susceptibil- 
ities of the impurity, which are different quantities that are 
positive-definite. 

In the high-temperature limit, the system behaves as the 
standard SIAM in the same limifc^: the impurity freely fluctu- 
ates between all four possible states, thus the entropy is In 4, 
while the magnetic moment is (0 + 1/4 + 1/4 + 0)/4 = 1 /8, 
since only two of the four states have magnetic moment of 
1/4. The charge moment TxQ,imp is 1/2 due to maximal 
charge fluctuations on the impurity. 

On the temperature scale T\ max(U, T), the directly 
coupled Majorana modes r\2a are frozen out (or, which is 
equivalent, the Dirac mode b is frozen out) and we enter a 
In 2 plateau in the impurity entropy which persists down to 
T = 0. The freezing-out of the charge degrees of freedom 
is also reflected in the effective charge moment TxQ,imp{T), 
which on the temperature scale T-y drops to some small (neg- 
ative) value and remains constant down to T = 0. The larger 
U is, the larger is the reduction of the charge fluctuations on 
the impurity, and the closer does the charge moment approach 
to zero. We emphasize that in the standard SIAM the effective 
charge susceptibility goes to zero at low temperatures and it is 
always positive. 

The residual In 2 entropy suggests that the other two Ma- 
jorana modes 7710- (i.e., the Dirac mode /) remain decou- 
pled; this is strictly true only for the non-interacting [7 = 
model, see Eq. ( fT2l . but an effective decoupling also occurs 
for U 7^ 0. The magnetic susceptibility curves indicate, how- 
ever, that the behavior within the In 2 entropy plateau is not 
trivial, see the two middle panels in Fig. [2] The most strik- 
ing feature is the strongly anisotropic behavior: while the ef- 
fective moments in the transverse x and z directions are re- 
duced ("screened") at low temperatures, the longitudinal (y- 
axis) effective moment reaches its asymptotic value on the 
temperature scale of T\ and remains constant down to T = 0. 
We also observe that the behavior is different depending on 
the value of the ratio U/T. (We remind the reader that this 
same ratio controls the stability of the non-magnetic solu- 
tion of the Hartree-Fock equations for the standard SIAM. 
For U/ttT < 1 the paramagnetic solution is stable, while 
for U /ttT > 1 the solution breaks the spin symmetry and 
the Hartree-Fock approximation is no longer applicable. This 
ratio thus determines whether the system is in the "Kondo 
regime" or not.) 

Let us first consider the case of large U, i.e., U/T > 1, 
which corresponds to the magnetic regime. The transverse 
magnetic susceptibility indicates a progressive reduction of 
the effective impurity magnetic moment. This reduction oc- 
curs at some low temperature scale T2 which appears to be 
exponential in U/T, similar to the exponential behavior of the 



Kondo temperature in the standard SIAM. The scale T2 may 
be defined, somewhat arbitrarily, through /x 2 .^) = ^ 2 (T = 
oo)/2 = 1/8. The temperature T2 is plotted as a function of 
U/T in Fig. [3] We find an excellent fit to the function 



T2 = c\U exp 



-C2 



(34) 



where c\ — 3.73 and C2 = 0.78 « 7r/4. This is dis- 
tinctly different from the expression for the Kondo temper- 
ature in the standard SIAM which takes the form Tk ~ 
U exp(— ttU/ST). The scale T2 appears to be in many re- 
spects similar to the low-energy scale studied in Ref. [3CJ, al- 
though the dependence on the parameters is exponential rather 
than algebraic. The physical mechanism which controls the 
scale T2 is not clear at present. 

40 1 1 
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Figure 3. (Color online) The temperature T2 where the effective 
transverse magnetic moment [1%. = Txxx,imp(T) is reduced to 1/8. 
The full curve is a fit which corresponds to Eq. 04b . 

Note that at zero temperature the impurity spin is not 
screened in the usual sense (which would be associated with 
a reduction of the impurity entropy to zero due to the forma- 
tion of a non-degenerate spin-singlet Kondo state). In partic- 
ular, the effective longitudinal moment is still non-zero, thus 
the spin is maximally anisotropic; it behaves as an Ising spin 
which cannot rotate away from the y direction. 

For small U of order T, we find somewhat different be- 
havior. In this case, the scales T2 and T\ are not well sepa- 
rated: the charge and spin susceptibilities are reduced simul- 
taneously. Again, the system has residual In 2 entropy and 
there is residual longitudinal magnetic moment. 

In analogy with the standard SIAM, we find that the interac- 
tion strength U plays the role of the effective bandwidth^, so 
that for constant U /T ratio the thermodynamic curves are sim- 
ply shifted to lower energy scales proportionally to U. Finally, 
it should also be remarked that at the particle-hole symmetric 
point, the low-temperature stable fixed point of the system is 
exactly the same irrespective of the values of the parameters 
U and T; again, a similar feature is also found in the standard 
SIAM^. 

It is interesting to note that the second scale T2 does not 
very clearly appear in the renormalization-group energy-level 
flow diagram, Fig. |H which is consistent with the constancy 
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Figure 4. Renormalization-group flow diagram. The NRG iteration 
step numbers corresponding to the temperature scales Ti and Ti are 
indicated. 



of Simp at low temperatures. The fact that the susceptibility 
curves nevertheless exhibit temperature dependence then im- 
plies that it is the nature of the low-energy excitations which 
changes with the temperature. There are no conserved quan- 
tum numbers in the problem, thus the information about the 
spin susceptibility is contained in the matrix elements, such 
as (i\S z \j), between the energy eigenstates. These matrix el- 
ements do have an RG flow. At the particle-hole symmetric 
point the system has double degeneracy (in fact the ground 
state and the low-energy excitations are all doubly degener- 
ate). For S = B = 0, this is a direct consequence of the 
system being non-ergodic, but the double degeneracy is also 
found for 5^0 (see also below), where the system is, in 
fact, ergodic. While the S = B = regime is in some sense 
pathological, we find that the behavior of the system varies 
smoothly as these parameters are varied from 0, therefore the 
results of the NRG calculations for S = B = are also phys- 
ically relevant in spite of the system being, strictly-speaking, 
non-ergodic. 



B. Thermodynamic properties: asymmetric case 

In the magnetic regime (U/T > 1), the departure from the 
particle-hole (p-h) symmetric point does not affect the behav- 
ior as long as 



\6\ < U/2 - cT, 



(35) 



where c is a number of order 1 , at which point the fixed point 
is destabilized, see Fig. [5] Notice that a similar condition 
was found in the bosonization approach [see Eq. d28l il even 
though the analysis was based on a large-r limit. This con- 
dition is similar to the region of existence of the "magnetic" 
local-moment fixed point in the standard SIA M 36 ' 39 . In the 
weak-interaction regime {U/T < 1), the magnetic fixed point 
is only stable in the immediate vicinity of the p-h symmetric 
point and the residual entropy is released at low temperatures 
even for small 5. There thus exists a two-dimensional sheet of 
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Figure 5. (Color online) Thermodynamic behavior away from the 
particle-hole symmetric point at 5 = 0. The arrows indicate the 
direction of increasing S. The transition point corresponds to |<5| ~ 
U/2 - cT with c of order 1. 



transition points in the (S, U, T) parameter space which sep- 
arates the regimes with or without the residual entropy. The 
quantum phase transition is of the Kosterlitz-Thouless type; 
the cross-over between the magnetic and non-magnetic fixed 
point occurs on a temperature scale which is an exponential 
function of 5 — S c , where 5 C is the critical value of the pa- 
rameter. A numerically determined phase diagram is shown 
in Fig. [6] We emphasize again that the existence of a phase 
transition is significantly different from the behavior of the 
standard SIAM, where the variation of the system properties 
in the (S, U, T) space is smooth and we merely move along a 
line of Fermi-liquid fixed points parameterized by the quasi- 
particle scattering phase shift. 

By changing the on-site energy 8, we tune the level occu- 
pancy n. The occupancy operator n is trivially related to the 
z-component of the isospin (also known as the axial charge) 
operator i z = l/2(n — 1), see Eq. 0. In this section, we have 
thus established that the magnetic fixed point is stable with 
respect to the non-zero isospin field i z . Due to the symme- 
try of the problem, this implies that the system is also stable 
with respect to moderately large isospin field i x , which can 
be induced by the superconducting proximity effect. We find, 
however, that non-zero isospin field i y will drive the system 
to the non-magnetic fixed point on the energy scale of i y . 
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Figure 6. Phase diagram in the (T, S) plane for fixed interaction pa- 
rameter U. 



C. Effects of the magnetic field 

In the presence of an external magnetic field, the entropy is 
released and the effective moment goes to zero at low temper- 
atures, see Fig. [7] The effect of the field is direction depen- 
dent. For a field along the "longitudinal" y direction, the en- 
tropy is always released on the temperature scale of B y . For 
a field along the "transverse" x and z directions the behav- 
ior depends on the U /Y ratio, in other words, it depends on 
whether the impurity is magnetic or not. For a magnetic im- 
purity (large U /Y ratio) the entropy is released on the temper- 
ature scale of \B±\, while for a non-magnetic impurity (small 
U/Y ratio) this only happens on a much reduced temperature 
scale, see Fig. [7] This behavior is described by the Kondo 
scale defined in Eq. (|29] l, although the parameters B± and 
r need to be rescaled by some constant factors in order to ob- 
tain full numerical agreement. 

The zero-temperature dynamical magnetic susceptibility 
curves for zero and non-zero magnetic field along transverse 
and longitudinal directions are shown in Fig. [8] In the absence 
of the field, the magnetic susceptibility in the transverse direc- 
tions diverges at lj = 0, while the longitudinal susceptibility 
goes to zero. All components of the magnetic susceptibility 
also feature a peak on the frequency scale of the atomic ener- 
gies (at ~ U, if we are in the Kondo regime, see Figs. [8]). In 
standard SIAM, the magnetic susceptibility peaks at Tr- and 
goes to zero at small frequencies as a linear function of lj, as 
mandated by the Korringa relation for Fermi-liquid systems^; 
the susceptibility curve also has a peak at lj ~ U/2. 

In the presence of a transverse magnetic field B x , all com- 
ponents of the magnetic susceptibility tensor peak on the scale 
of the magnetic field and then go to zero. This is even true 
for the longitudinal susceptibility which in the absence of the 
field goes to zero. For longitudinal magnetic field B y , the 
transverse components of the magnetic susceptibility peak on 
the scale of the magnetic field, while the longitudinal suscep- 
tibility peaks on the same scale as in the absence of the field 
and no further peak emerges at lj ~ B y . These features of 
the magnetic response of the system are characteristic for im- 
purities coupled to helical Majorana edge states and clearly 
distinct from that of the standard SIAM. We also note that the 
results do not seem to agree with those found for an equiva- 
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Figure 7. (Color online) An impurity in an external magnetic field 
along x, y, and z directions, respectively, with constant strength |B|. 
For each direction of the field, the parameter U is swept from the 
weak-interaction to the strong-interaction (Kondo) regime. The ar- 
rows indicate the direction of increasing U. 



lent Kondo impurity problem in Ref. |27|. The difference can 
be explained in part by the fact that the Anderson and Kondo 
impurity models are not fully equivalent; see Sec.lVl 



D. Spectral functions 

The zero-temperature spectral functions of the impurity are 
shown in Fig. [9] In addition to the charge-fluctuation peaks at 
lj ~ U/2, as in the standard SIAM, one observes a sharp res- 
onance on the scale of with an inverse-power-law shape, 
i.e., the spectral function diverges. (Of course, one cannot 
attach a scale to a power-law function. Instead, T2 roughly 
corresponds to the energy where the cross-over to the power- 
law behavior occurs.) This inverse-power-law divergence re- 
places the Kondo resonance of the standard SIAM, which at 
the lowest energy scales looks like a parabolic peak and is fi- 
nite. For strictly decoupled Majorana modes one would ex- 
pect a delta peak at zero frequency; the "broadening" into 
the inverse-power-law peak is thus an interaction effect re- 
lated to the Anderson orthogonality catastrophe physics. The 
exponent depends on the interaction strength, see the inset 
in Fig. [9] We find that the exponent a is well fitted by 
a = 1 - I arctan((7/4r). 

We also observe that the anomalous spectral functions are 



9 



Transverse field 



— B =10" 



B =10" 




Figure 8. (Color online) Dynamical magnetic susceptibility func- 
tions for an external magnetic field applied along the transverse (x) 
direction (top panels) and along the longitudinal (y) direction (bot- 
tom panels). 



non-zero, which is expected for an impurity coupled to a su- 
perconductor. 
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Figure 9. (Color online) Spectral properties of the impurity. The 
non-zero spectral functions are shown. We introduce the notation 
do = d\, d\ — dt, d-2 — di, ds = d^ and define the spectral func- 
tions as Aij(uj) = —(l/n)Im((di;dj)) UJ +is- Non-zero anomalous 
spectral functions A02 and A13 thus indicate the presence of pairing 
correlations. The inset in the bottom panel represents the power-law 
exponent as a function of the interaction strength. 



V. SCHRIEFFER- WOLFF TRANSFORMATION 

We now consider a generalization of the Hamiltonian H de- 
fined in Eqs. ([T]|5]l. We take the conduction-band Hamiltonian 
as in the standard SIAM: 



H 3 = X! efeC L c ^- 



(36) 



This conduction band has twice the number of the degrees of 
freedom as that described by Eq. © since there are now two 
Majorana modes at each k: r/kia and rjk2a- The Majorana and 
Dirac modes are related by 



T}k\o 



1 / t 
^= [Vkla + lVk2a) , ko 

1 

= 7i 



Ckc 



-k<j ) ' 



1 
1 



klc 



Cfe CT 



irik2a) 



(37) 
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The hybridization term allows for unequal coupling of the Ma- 
jorana components: 



H' 2 = a^V k (ri kla d a + 4»/fcla)/v^ 
+ ft Y, Vk{~irjk2<yd a + id\r\u2cr) I V2, 



(38) 



where a and ft are some real numbers. For a = ft = 1, this 
is the standard SIAM, while for a = 0, ft = 1 we recover the 
Majorana SIAM studied in this work (up to some differences 
in the notation and a factor of y/2 in the definition of V k ). We 
can thus smoothly interpolate between the two limiting cases. 
In the following, for reasons of simplicity we assume Vk to be 
independent of k, i.e., V k = V . 

We now proceed to perform the Schrieffer- Wolff 
transformation^^ for the generalized Anderson impu- 
rity model in the large-C/ limit in order to derive an effective 
Kondo model. The low-energy subspace consists of those 
many-particle states where the impurity is singly occupied, 
and we denote by Pq the projection operator onto this sub- 
space. Then Qq = 1 — Pq is the projector onto the orthogonal 
("high-energy") subspace. We introduce the superoperator 



O(X) = P Q XQ + Q XP . 



(39) 



We find that 0(H' 2 ) = H' 2 , thus H' 2 is a block-off-diagonal 
operator. We next introduce the superoperator 



(jlgPOjj) 

E; — Ei 



\i)(j\-n.c. 



(40) 



Here {\i}} is an orthonormal basis of Hq = H± + such 
that H \i) = Ei\i). We assume that the states \i) in Eq. d40l 
belong to the low-energy subspace, while \j) belong to the 
high-energy subspace (the — H.c. term then generates the re- 
maining terms). To lowest order in hopping aV, one has 



and 



We find 



s = c{H' 2 ) + o(y 2 i 



H eS = H P a + -P [S, H 2 ]P + 0(V 4 ) 



(41) 



(42) 



C{H' 2 ) = Y 

- • }•{'■, } 

k{*i}>(7{*,-}| 



(a{fc i; }|g 2 | 7 {fc 3 }) 

E <r{ki} ~ E y{k 3 } 
- H.C., 



(43) 



where a € {f, 4-} indexes the impurity states in the low- 
energy subspace, 7 6 {0, 2} indexes the impurity states in 
the high-energy subspace, while {fcj} and {kj} are the occu- 
pancies of the conduction-band states. All terms in H 2 are 
such that the sets {ki} and {kj} must differ by the occupancy 



of a single level. We obtain 

S = ^^^(l-n_ CT )cL4-H.c. 

fc<7 

ho 
ho 

+£ 



{ft 


- a)V/2 




£fe + e 


(a 


+ ft)V/2 




(k - e 


W 


- a)V72 




+ £ + £/ 


(q 


+ 0)V/2 


£k 


-e-U 



H.c. 



(44) 



n~aC k(J do 



H.c. 



The effective Hamiltonian is very complicated. A number 
of terms are similar to those in the Schrieffer- Wolff transfor- 
mation for the standard SIAM, Ref. |4lj, but with coefficients 
which depend on a and ft (and reproduce the standard val- 
ues in the a = ft = 1 limit). Furthermore, some magnetic 
anisotropy also arises. For example, the exchange-coupling 
terms can be written as 

#ex= Y Y^k'^^d^d), (45) 

i£{x,y,z} kk' 

where S = {S x , S y , S z } = (l/2)<r are the Pauli matrices, 



_ /Cfcf 



CfeJ, 



and *&, 



(46) 



are field operators, and the effective Kondo exchange coupling 
coefficients Ji are equal to 



J T. J Z 



V 2 (a + ft) 2 /2 V 2 (a-ft) 2 /2 



r2f„, 1 o\2 



(47) 



J„ = 



_ F 2 (a + /j)72 _ V 2 {a-ftY/2 

e k -e-U e k + e + U 

V 2 (a + ft) 2 _ V 2 {a-ft) 2 

e k - e e fc + e 

_ y 2 (a + ^) 2 F 2 (a-/3) 2 
e k -e-U e k + e + U ' 



Only in the a = /3 = 1 limit is the effective exchange scatter- 
ing isotropic. At the particle-hole symmetric point we find 



J x J z 



Ju 



2V 2 



2V U 



a \e k + U/2 e k -U/2 / 
(a 2 +ft 2 ) ( 2V 2 2V 2 



2 \e k + U/2 e k -U/2 
The direct s-d interaction is 

J: 

W - 

kk 

where 

V 2 (a + ft) 2 /A V 2 {a-ft) 2 /4 



(48) 



W 



£fc - e 



(49) 



(50) 
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The term which can be absorbed in the Hamiltonian Hq 
Hi + H' s can be written as 



ha 



Ho = -J2( W ' + l J ' n d,-<r)n d , 



(51) 



where 



IT" 



J' 



, V 2 (a + /3) 2 /4 V 2 (a-(3) 2 /2 



£fe - e (k + e 
V 2 (a-/3) 2 /4 

e k + e + U 

V 2 (a + f3) 2 /2 V 2 (a-pf 



(52) 



£fc - e (k + e 

V 2 (a + f3) 2 V 2 (a-(3) 2 



e k - e - U e k + eu 
The two-particle hopping term is 

H ch = - ^2 hhc)., _ a c\. a d rT d- a + H.c, (53) 



kk f a 



with 



'2h 



V 2 {a + p) 2 V 2 {a + (if 



ek - e 



ek - e-U 



(54) 



These are, however, not all the terms which appear in the 
effective Hamiltonian. Rather than enumerate all remaining 
contributions, we restrict our attention only to those which 
maintain the single occupancy of the impurity orbital. They 
can be written as: 

H *i = L \ c U c i'l ~ c k'i c kf) ( d \ d ± ~ d l d V ( 55 ) 



fcfc' 



with 



V 2 (P 2 -a 2 )/8 V 2 (a 2 -f3 2 )/8 



ek - e 



ek + e 



+ 



V 2 {p 2 -a 2 )/8 V 2 ((3 2 -a 2 )/8 



(56) 



e k - e - U 



e k + e + U 



We now discuss the a = 0, ft = 1 limit, which is relevant 
for discussing an Anderson impurity coupled to the edge states 
of a helical topological superconductor. The most important 
terms in the effective Hamiltonian are those which affect the 
impurity spin degrees of freedom: 

H c g = H cx + H si . (57) 

The exchange coupling constants are anisotropic: 

V 2 /2 V 2 /2 V 2 /2 V 2 /2 



J 



e k - e ' e k + e ek — e — U ek + e + U' 
V 2 /2 V 2 /2 V 2 /2 V 2 /2 



J' II 



ek - e ek + e ek - e-U ek + e + U 



(58) 



At the particle-hole symmetric point this simplifies to 



J-u 



0, 

AV 2 U 
U 2 - 4e k 



W 2 

IT' 



(59) 



where we have used ek — > 0. Within the same approximation, 
we also have 



L 



v_ 

u 



We combine all terms and write 



4:SySy 



Yl 
u 

w 2 

jj \Sy — Iy) s y 



(60) 



(61) 



Here upper-case operators / and S correspond to the 
conduction-band isospin and spin, while the lower-case op- 
erators are those of the impurity. Furthermore, I + = I x + il y , 
etc. We thus conclude that in the large-C/ limit, the impurity 
spin degree of freedom couples to a mixed S y — I y mode, 
which can be expressed in terms of a single Majorana chan- 
nel. This result was postulated without derivation in Ref. [27l 
It is important to note, however, that the derivation and the 
effective model only make sense in the large-?/ limit, i.e., the 
coupling constant Jk — W 1 jU verifies Jk <C V <C U. 
The phase transition discussed in Ref. |27| occurs at a value of 
pjfc ~ 0(1) which is therefore unphysical. Thus we con- 
clude that such phase transition is not expected in real mag- 
netic impurities coupled to topological superconductors. 



VI. CONCLUSION 

We have studied the properties of the modified single- 
impurity Anderson model where the impurity couples only to 
half of the degrees of freedom of standard fermionic particles 
in the conduction-band continuum, i.e., to Majorana fermion 
channels. We have shown that the model is related to the two- 
level resonant model with attractive charge-charge interaction 
and to the (antiferromagnetic or ferromagnetic) anisotropic 
Kondo model. Two different stable low-temperature fixed 
points have been identified: one corresponds to magnetic im- 
purities and is characterized by residual magnetic moment 
and entropy, the other corresponds to non-magnetic impurities 
with no residual degrees of freedom. The phase diagram sep- 
arating these two regimes has been established. We have also 
shown that the magnetic field always destabilizes the magnetic 
fixed point, however the response is strongly anisotropic. 

For magnetic impurities, the residual degree of freedom 
corresponds to the operator s y + i y . While the charge de- 
grees of freedom are quenched due to electron-electron re- 
pulsion, for impurities which are not in the extreme Kondo 
limit (U 3> T) the isospin degree of freedom determines the 
quantitative aspects of the problem. For this reason, we have 
shown that it is important to properly map the original Ander- 
son model to an effective Kondo model, rather than start by 
postulating a Kondo-like model. 
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Finally, we would like to stress that the predictions we have 
made can be put to an experimental test since dynamical mag- 
netic susceptibility is in principle measurable experimentally. 



and obtain 



H 2 = J2\ V ^ ir 1^2a- 



(A4) 
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Thus we find, again, that only two Majorana local modes ^2<j 
are coupled to the impurity, while the remaining two modes 
£icr are decoupled. If £i^£u i s rewritten in terms of the orig- 
inal Dirac operators, it is found that it corresponds to a linear 
combination of spin and isospin operators in the (xy) plane: 



Appendix A: Complex hybridization matrix elements 

For completeness, we now study the case of general com- 
plex and spin-dependent hybridization matrix elements, that 
is, the coupling Hamiltonian is written as 

H 2 = Y. {y^kcda + VZdUkv) , (Al) 
fcfj 

where K- is a complex number, V a = \ V a \ cxp(i9 a ). 
Introducing r\i a as in Eq. ©, we obtain 

H 2 = ^2 \V a \V2i (cosd a ri ka T]2a + smO^k^^) . (A2) 

ka 

We now perform a change of basis to 



»Cl4.flt = ^(fy - 8t) S X + COS(6» 4 - 9^Sy 



= cos 6*0.771,7 - sin9 a j]2a, 
£ 2 <r = sin^ryio- + cos9 a rj 2c r, 



(A3) 



(A5) 



- sin(0 t + 6±)i x + cos(6» t + 9^)i y . 
In the absence of spin-dependence, i.e., for 6f = 6±, we obtain 

*6u.fit = s y _ sm(26)i x + cos(29)i v , (A6) 

which implies that the privileged "longitudinal" magnetic axis 
is still along the y direction, however the privileged isospin 
axis deviates by an angle of 29 from the y-axis. For 9^ — 
&l = we recover the result from the main text, 



*£u£it = Sy+i y . 



(Al) 



We remind the reader that the z-axis is defined by the spin 
polarization of the Majorana edge modes. The x and y-axes 
are thus always defined relative to this spin quantization axis. 
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